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Treating problems in full general relativity is highly com- 
plex and frequently approximate methods are employed to 
simplify the solution. We present comparative solutions of 
a infinitesimally thin relativistic, stationary, rigidly rotating 
disk obtained using the full equations and the approximate 
approach suggested by Wilson & Mathews. We find that the 
Wilson-Mathews method has about the same accuracy as the 
first post-Newtonian approximation. 
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I. INTRODUCTION 

The numerical solution of Einstein's field equations in 
three spatial dimensions in simulating for example the 
merging process of two neutron stars is very complicated, 
and only very recently have methods been developed for 
treating such problems in full general relativity. One 
of the main difficulties is encountered in the treatment 
of the radiation degrees of freedom. An approximation 
scheme developed by Wilson & Mathews JO simplifies 
the problem by imposing a conformally flat condition on 
the three metric, throwing away all radiation degrees of 
freedom. This reduces firstly the number of field equa- 
tions and secondly reduces them to Poisson like equations 
exclusively. In general, the first post-Newtonian approx- 
imation of Einstein's theory is identical to the one of 
the Wilson-Mathews scheme, and in spherical symmetry 
Wilson-Mathews is identical with Einstein. 

As pointed out by Cook et al. |Q, this approximation 
should be viable and verifiable for systems in a true equi- 
librium state. They investigated the accuracy of Wilson 
& Mathews' approximative scheme on isolated rapidly 
rotating relativistic stars, and found largest errors of 
about 5% for any quantity analyzed. As the approxi- 
mate scheme has been used in calculating the evolution 
of a binary system || which is very far from spherical 
symmetry, it is desirable to have test cases, that depart 
stronger from spherical symmetry than rotating stars do. 

Here we present the model of an infinitesimally thin, 
pressureless disk that rotates rigidly. This problem was 
considered by Bardeen & Wagoner Q up to the tenth 
post-Newtonian approximation and has been solved ex- 
actly by Neugebauer & Meinel S. Thus, it serves as 
an ideal test bed to check the validity of any numerical 
method solving the stationary, axisymmetric field equa- 



tions H and in particular Wilson & Mathews' approxi- 
mation scheme. To demonstrate the general accuracy of 
the numerical method employed we present the numerical 
solution of the full Einstein equations as well. 



II. BASIC EQUATIONS 

To describe the stationary, axisymmetric equilibrium 
of an infinitesimally thin disk we use a cylindrical coor- 
dinate system (c = 1, c velocity of light) 



X{) = t, X\ = p. X2 = Z, X3 = (f 



(1) 



where the disk is located in the equatorial plane (z=0), 
and the rotation axis coincides with the z-axis. 



A. Full Einstein equations 

The general line element for an axisymmetric and sta- 
tionary mass configuration can be written as Q] 

ds 2 = e 2 » (dp 2 + dz 2 ) + \pBe-"{dtp - todt)} 2 - e 2u dt 2 , 

(2) 

where the metric potentials v,h,uj and B are functions 
of p and z. As the disk has no vertical extension, the 
global solution of the problem may be obtained by solving 
the vacuum field equations, where the presence of the 
disk manifests itself through jump conditions in the first 
(vertical) z— derivatives of the metric functions. 

As shown by Q and 01, the vacuum Einstein equations 
Gij = for the metric (|2|) are then given explicitly by 

V(Vf) = -p 2 B 2 e^ u Vio ■ Vlu - —VB ■ Vv (3) 
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(6) 



with the V-Operator in cylindrical coordinates (p, z), and 
the Cartesian (in p and z) Laplace operator A c . 



f 



The boundary conditions in the equatorial plane within 
the disk region p < pd, where pd denotes the outer disk 
radius are then given by (see H) 



Here a and P denote the surface density and pressure, 



and v v = pBe 2u (fl 



is the orbital velocity of the 



gas. We note that in case of a pressureless (P = 0) dust 
disk the metric function B is identical to unity Q. 

Radial hydrostatics yields the relation e v (l — vVj ^ 2 = 
const, in the disk, where p < pa- In the outer part of 
the equatorial plane, p > pd, and along the polar axis 
symmetry conditions are used for the metric coefficients. 
At infinity p, and u tend to zero, and B = 1. 



For our test case of a rigidly rotating dust disk (B = 1) 
we may summarize the exact Einstein metric and the 
approximative Wilson-Mathews form as 



(17) 



(18) 

respectively. Comparing ( |l7| ) and ([l^), it is striking 
that the only difference between exact and approximate 
method for the dust disk example lies in the function 



v tZ = 2na f 

1 - vfp 


(7) 


ds 2 = e 2 » 


(dp 2 + dz 2 ) 


^p 2 e- 2 ^(d^-Lodt) 2 - 
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and 






B jZ =0 

P,z = — 27TCT. 


(9) 
(10) 


ds 2 = e 2 * 


[(dp 2 + dz 2 ) 


+ p 2 (d V -Codt) 2 } -e 2i/ dt 2 , 



The function £ gives the deviation from 
three-dimensional conformal flatness in Einstein's theory 
and it has been shown that in the first post-Newtonian 
approximation £ vanishes for the dust disk solution H . 



C. Method of solution 



B. Wilson-Mathews formulation 

The method developed by Wilson and Mathews (j]] 
consists in approximating the line element (^) by using a 
spatially isotropic form and neglecting the non-isotropic 
parts. Specifically, for the problem at hand we rewrite 
the exact metric as 

ds 2 = e 2fi [(dp 2 + dz 2 ) + p 2 (dip - Codt) 2 ] 

+ p 2 B 2 (dip-Codt) 2 -e 20 dt 2 (11) 

with the functions B, p, and Co. 

By comparison with (0) we may identify 



(12) 



and for the metric function B we obtain 

B 2 = (e 2 ^ + B 2 )e 2i> . (13) 

Applying the approximation by Wilson & Mathews, we 
neglect the non-isotropic part and set B = 0. Simulta- 
neously we have to omit the corresponding field equation 
(for B) 7 and we obtain the approximate field equations 
from (|), §), and (§) by setting B = 

V(V£) = ^p 2 e 2 ^-^ Vw • VQ - V(£ + p)-Vv (14) 



V(p 2 Vuj) = -p 2 \7(3p - v) ■ Vu>, 



(15) 



1 i C\~ 

A c p = -p 2 e 2 ^-^VCj ■ VCj + Vp ■ VP + (16) 
4 pop 

Furthermore, we have to substitute for B in the equation 
for the orbital velocity, and obtain v v — pe^~ v \£l — Co). 



The field equations for the exact and the approximate 
case together with the boundary conditions are solved 
numerically for the problem of a rigidly rotating dust 
disk. As the numerical procedure has been described in 
more detail by Kley || we give here only a brief sketch 
of the method. To prevent problems with the conditions 
at infinity the computational domain [0 : oo; : oo] is 
mapped onto the unit square where the radius pa of the 
disk is located at 1/2. This domain is covered here by 
50 x 50 gridcells and the partial differential equations are 
discretized to second order spatial accuracy on this grid 
yielding a coupled system of equations. Together with 
the boundary jump conditions these equations are then 
solved iteratively until convergence has been achieved. 
The accuracy of the solutions may be estimated using 
the analytical solution J5| and post-Newtonian approxi- 
mation fE| of the dust disk. 



III. RESULTS 

To give an idea of the degree of non-isotropy of the 
three metric for the exact Einsteinian disk we display 
in Figure 1 the behaviour of the function £(p) within 
the disk for different values of the central redshift. In 
the figure the numerically obtained solution of the field 
equations (|))-(^|) together with the boundary conditions 
(f?)-(|l0|) is compared to the exact solution as given by [||. 
The errors over the whole disk for all central redshifts 
are usually well below 1% and only at the outer edge 
of the disk the worst deviation is about 2% for z = 3. 
This demonstrates the accuracy of the numerical method 
employed, even for the moderate gridsize used. 
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FIG. 1. The function e 2C which gives the deviation from 
local flatness plotted for the Einsteinian dust disk as func- 
tion of the normalized radius p/pd for three different central 
redshifts of the disk. The diamonds refer to the numerical 
solution of the full Einstein equations, while the solid lines 
refer to the exact solution. 
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FIG. 2. The metric function u as a function of the nor- 
malized radius p/pd for three different central redshifts of the 
disk. The diamonds refer to the numerical solution of the full 
Einstein equations, while the crosses denote the numerical so- 
lution obtained by using the Wilson-Mathews approach. The 
solid lines refer to the exact solution. 

We compared the numerically obtained approximate 
metric coefficients v 1 fl and Co versus radius with the exact 
solution and find that in case of the "lapse" potential v 
the deviation from the exact case is always below 2.5%, 
even for central redshifts as large as 3. For the function 



H the deviation is larger, 2.3% for z = 0.5, 5.5% for 
z = 1.0, and about 20% for z = 3. The results for the 
metric function u> are displayed in Fig. 2. The errors are 
very large for higher redshifts, about 9% for z = 1 and 
up to 48% for z — 3. The approximate solution (crosses) 
shows also the wrong trend of the evolution of u) with 
redshift. While in the Einsteinian disk case u> increases 
first with redshift, and then decreases for z > 2, in the 
approximative solution u> seems to increase steadily with 
z. 
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FIG. 3. The quantity QA/o as a function of central red- 
shift z for the dust disk. The open squares correspond to the 
numerical solution of the full Einstein equations, the crosses 
to the approximate Wilson-Mathews procedure. The solid 
line refers to the exact solution, and the dotted and dashed 
line to the first and second post-Newtonian approximation 
respectively, of the Einstein equations. 



The overall deviations of the metric functions from 
the exact case may be understood qualitatively by look- 
ing at the line elements (17) and (18|). As the lapse is 
not too much affected by the simplification in the spa- 
tial three metric it will be calculated quite accurately 
in the Wilson-Mathews approach. The rotational coef- 
ficient uj is affected mostly and hence shows the largest 
errors. We may nevertheless infer that at least for central 
redshifts up to z — 0.5 the approximate solution is very 
accurate, but for higher relativistic disks the agreement 
deteriorates. Note, that z = .55 was the maximum (po- 
lar) redshift tested in the comparison for rotating stars 
by Cook et al. @. 

As a further illustration of the range of applicabil- 
ity of the Wilson-Mathews approximation, it is interest- 
ing to plot global quantities. In Fig. 3 we display the 
dependence of the dimensionless quantity QMq versus 
the parameter 7 = z/(l + z). One notices again the 
very good agreement of the numerical and exact solution 
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of the full Einstein equations. Interestingly, the post- 
Newtonian approximations approach the exact solution 
from below while the Wilson-Mathews approximation lies 
always above the exact solution. In case of the full Ein- 
steinian disk for infinite central redshift £1 * Mo reaches 
a finite maximum of about 0.8. 
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FIG. 4. The binding energy as a function of flMo. The 
symbols are identical to Figure 3. 
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FIG. 5. Dimensionless combination of the parameters 
measurable at infinity: gravitational mass M, angular mo- 
mentum J and angular velocity £1. The symbols are identical 
to Figure 3. 

In Fig. 4 the normalized binding energy E^/Mq = 
1 — M/Mq is plotted versus SIMo. For a given rest mass 
Mo and angular rotation rate f2 the Wilson-Mathews ap- 
proximation yields a lower binding energy for the dust 



disk than the Einstein disk. The dotted line for the post- 
Newtonian solution ends at O * Mo = 0.57 because this 
refers already to z = oo (cf. Fig. 3). 

Finally, in Figure 5 the dimensionless rotation rate MO 
is plotted versus the dimensionless total angular momen- 
tum M 2 /J, parameters that are measurable at infinity. 
The value 20M = M 2 /J = 1 corresponds to the ex- 
treme Kerr limit of the dust disk solution (see ||). Again, 
the Wilson-Mathews approximation lies beyond the post- 
Newtonian approximation and does not yield the correct 
black hole limit. 



IV. CONCLUSIONS 

We have calculated numerically model sequences of a 
rigidly rotating disk of dust using a) the full Einstein 
equations and b) the modified Wilson-Mathews scheme. 
By comparing those to the known exact solution by 
Neugebauer & Meinel (j| , we have shown that the Wilson- 
Mathews approximation yields results which are accurate 
to about 5% if the central redshift of the disk is smaller 
than 0.5; but may be significantly higher for stronger rel- 
ativistic models. We conclude that the Wilson-Mathews 
approach for solving the Einstein equations yields results 
which are at most as accurate as the first post-Newtonian 
approximation, at least for rigidly rotating disks of dust. 
Hence, its application to stronger relativistic configura- 
tions e.g. in close, compact binary systems has to be 
handled with caution. 
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